X-602- 74-166 

PREPRINT 


r* NASA TM X- yo 
THE AVERAGED MOTION 
OF A CHARGED PARTICLE 
IN A DIPOLE FIELD 


S ' 

A. J. CHEN 
DAVID P. STERN 


(NASA-TM-X- 70668) THE AVERAGE MOTION OF 
A CHARGED PAPTICXE IN A DIPOLE FIELD 
(NASA) 15 p HC $4.00 CSCL 2QH 

G3/24 

t 

MAY 1974 


■ GODDARD SPACE FLIGHT CENTER 

GREENBELT, MARYLAND 



874-26185^ 

Unclas 
40939 j 


For information concerning availability 
of this document contact: 

Technical information Division, Code 250 
Goddard Space Flight Center 
Greenbelt, Maryland 20771 


(Telephone 30b982-4488> 



'I - 


The A v 
P a 


e r a g e d 
r t i c 1 e 


Mot ion 
in a D 


of a C 

i p o 1 e F 


h a r g e d 
i e 1 a 


A. J. Chen* and David P. Stern 

Goddard Space Flight Center 
Greenhelt, Maryland 20771 


* HAS/lffiC Research Associate 



2 


INTR ODUCT I OK AND NOTATION 


The average motion of a charged particle trapped in the geomagnetic 
field can usually he represented with adequate accuracy hy assuming the 
conservation of the first two adiabatic invariants £e*g. Northrop , 1265] 


/ 


p = P f / 2 m B (1) 

J p H ds (2) 

where integration is along a field line between mirror points „ In addi- 
tion, specification of the averaged motion also involves the parameters 
defining the magnetic field line to which the particle is attached. This 
role is usually served by the Euler potentials ( ot , P> [ e.g. Stern , 

197 Q ] ; given (j, |_l) of the particle and ( , P } of the field 

line to which it is attached at some initial time t^ , the "drift surface" 
described by Its averaged motion is completely determined. 

In the absence of electric fields and time dependence the energy is 
conserved and the drift surface may be described by properties of the 
trajectory not involving the energy, such as the mirroring field B m and 
the quantity I' = j/v . In this case and the initial values of 

) are sufficient to define the drift surface, since I' may be 
derived from them. 

If in addition the magnetic field has the dipole configuration it 
is often convenient to use instead of the sine y of the equatorial 

pitch angle. An appropriate set of Euler potentials is then, in polar 
coordinates 
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= ( a 2 g° / r ) sin 2 © (3) 

l 

P = a'f (4) 

where a is the earth's radius and g° = -0.31 gauss is the dipole 

l 

harmonic coefficient „ In many applications it is convenient to introduce 
a related quantity 

L = a g° / o4 (5 ) 

i 

equal to the equatorial distance , in earth radii, of the field line 
with the given value of © L (one could use L as an Euler potential, 
hut the form of & then is not as simple). With the above definitions one 
can define the dimensionless "longitudinal integral" 


B 

m 



B' 

m 


evaluated between mirror points and proportional to I' . Here the initial 
values of © L and of y are all that is needed to define the drift surface, 
since I depends on y and is absent due to axial symmetry. Appro- 

ximations of the function l(y) have been derived by- Schulz [Mi] and 
by others* 


In the geomagnetic field this approach is often use$ at high energies , 

where the effects of electric fields can he neglected© However, at 

energies of the order of one kilovolt this method no longer works, since 

t&^lnergy may change appreciably : in this range \i and J , rather than 

B and I 1 (or l) are to he regarded as conserved quantities* The pur- 
m 

pose of this work is to provide the means for treating this more general 
motion, assuming a dipole field© 
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THE AVERAGED HAMILTONIAN 


The chief tool in this treatment is the averaged Hamiltonian K(<J, p,«C, P / 
the properties of which have been developed "by Northrop and Teller t 1260 , 
eqs. 30 ; see also Northrop , 1963 . eqs. 3.68 ] . Assuming the adiabatic 

invariants J and p are conserved and denoting by W the kinetic energy 
(in all what follows the nonrelativist ic form m v 2 / 2 will be used), the 
averaged Hamiltonian has the form 

K - W(«t, J, (i) + e^ + (e/c)ct ^P/'St (7) 

where W is independent of P due to the axial symmetry of the dipole 
field. For adiabatic ity to hold, the electric field must be orthogonal 
to B at least to the lowest order; if such orthogonality is assumed 
and no time dependence exists, the electric potential $ will be a 
function of (o4,P ) only and in general, the part of K not including 
W will depend on oL and p only [~ Northrop and Teller , i960, eq. 9 J . 

Two limiting cases are easily calculated. If J = 0 the particle is 
confined to the equatorial plane and one gets 

W ( oi , 0, |i) = p B = p jg^j/ L 3 - P(p, L) (8) 

This form has been widely used in examining the motion of charged parti- 
cles in the equatorial plane of model magnetospheres T Chen , 1970 ; Stern , 


On the other hand, if p = 0 then v f/ = v . In that case , if 
S denotes the total length of the guiding field line, equation (2) 
reduces to 

J = 2mvS = 2 m v 1 ^ a L ( 9 ) 
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where 

I 0 = 2 + 3“ ^ ln( 2 + 3^ ) = 2.76 

Expressing the kinetic energy gives 


W (ct, J, 0) = (j®/8 m Ij a 2 ) L" 2 = Q(J, L) (10) 

Here an approximation will he derived for W which can he used over 
the entire range of parameters. In doing so, account should he taken of 
the special nature of the dipole field, which constrains the allowed form 
of W . 

The dipole field has the property that all it-s field lines have the 
same configuration - only the scale length of the configuration varies* 

This reduces the functional form of W from a dependence on 3 variables 
to a dependence on only two* 

Por proof, let y he introduced as auxiliary variable* Then from (l), 
(6) and (10) 

r 2 (y) = Q ij / w (11) 

y 2 = P / W ( 12 ) 

If y is eliminated between (ll) and (12) a relation is obtained of 
the form 

P / ¥ = F(Q / I 2 W) (13) 

If W is extracted from this, it can depend on p , J and oL only 
through the functions P and Q ♦ A convenient representation of it, with 
the correct dimensionality and limits, would he 

W(P, Q) = P + Q 4- \ pSk Q 1 " Sk U4) 

where 0 < 1 , since Q vanishes when W = P and vice versa. 
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It simplifies the calculation if Q is replaced by 

X = ( Q I* / P ) ^ - I (y)/y ( 15 ) 

This reduces (14) to the form 

W(p, X) = P ( 1 + X 2 /^ + X % ) (16) 

Comparison with (12) shows that the expression in parentheses simply equals 

_ 0 

y . The results derived here will he presented in this notation. 

THE APPROXIMATION TO W 


One term in the expansion, the term linear in X can he derived 
from the harmonic-oscillator approximation developed hy Schulz [12a] 
for values of y near unity 

l(y) = 0.7405 ( 1 - y 2 ) (17) 


If y = 1 - $ , then 


X/O.7405 = 


Substituting in (12 ) gives the coefficient of X near y = 1 as 1/0.7405 
(the term proportional to X 2 may he neglected there ) 

W = P/ y 2 = P ( 1 + 2 S ) 


p ( 1 -t- 1.3505 X ) 


( 18 ) 
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A relatively accurate approximation to W , including two terms in 
addition to the three developed here, is 

w = p ( 1+ 1.3505 x - 0.030425 x 1 * - / 5 * 0.10066 x 5//5 + X 2 / 1^ ) 

= P f(X) ^ 1 9) 

As a test ve note that by (12) this should satisfy 

y 2 f(x) - 1 (2°) 

Actual values of this product, as functions of the mirror colatitude 

ft , are listed in Table 1 . The largest errors fall in the range which 
m ’ 

in most practical cases is inside the loss cone and for most of the iield 
line the error is 1%> or less* 

APPLICATION 

The relation (20) establishes a certain relation betveen X and y , 
or equivalently betveen X and the equatorial pitch angle ^ ( y - s ^- n ^ 

and this relation is plotted in Figure 1 « This plot can be used to deter- 
mine the energization of particles drifting earthvard in a combined 
electric and magnetic field, in the folloving vay o 

Suppose a particle starts at L = Li with an equatorial pitch angle 
and arrives at L = L2 • ve wish to know by how much has its kinetic 
energy changed „ The initial kinetic energy is readily derived from the 
given data, for by (12) 


W 


\x g° j L 3 sin 2 A 


( 21 ) 
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From Figure 1 the initial value Xf corresponding to is read 

off. By (15), ( 8 ) and (10) 

X = (j^/a 2 8 m(j.|g^| )s L 2 = constant • L 2 (22) 

Therefore the final value of X is 

Xa - Xj. (WLi) a (23) 

With X 2 having "been derived, the corresponding value of the pitch 
angle is readily read off Figure 1 , after -which (21 ) is used 

to obtain the final energy W 2 . 

From ( 8 ) and (10) it is evident that in the limiting cases 

-3 -2 

J = 0 and |i = 0 W increases as L and L , respectively, and 
near these limits the above more elaborate method need not be used. It 
becomes useful in the intermediate range where neither of the adiabatic 
invariants can be considered to be small . 

MEAN DRIFT 


By the formula developed by Northro p and Teller | ^1960 j the average 
motion of the particle satisfies 

<o4> = -(c/e) ^ K/^f 3 (24) 

* 

= (c/e) W'K ( 25 ) 

The differentiation of the electromagnetic terms in equation ( 7 ) is 
straightforward and therefore only that of W will be discussed. In the 
special case of the dipole magnetic field no dependence on p exists, so 
there is no contribution to < U. > from this term. For deriving ( 25 ) 
must substitute 


one 



- 9 - 


1 .2 

X = J (8 .a m vx) 2 oi 2 

(26) 

P = li ci 3 / a 3 (g °) 2 

l 

after which differentiation is straightforward „ In order to give the 
motion an intuitive representation it is useful to map it into the equatorial 
plane, where every point corresponds to a pair of values of (o£, P ) 
and can therefore represent the field line passing through it. Let V he 
the velocity of the point at which the field line to which the particle 
is attached crosses that plane .Then if — 0 and all variables are 

evaluated in the equatorial plane 


<<*> = V-V* 

= v-vp 


(27) 


Substituting (27) in ( 24 -) and (25) and subtracting the former from 
the latter then gives 


VK = (c/e) V X B 


(28) 


Since V is in the equatorial plane, it has no component parallel to B , 
from which follows (compare eq. 51 , Morthro p and Teller , i960 ) 

V = (e/cB 2 ) B X 7 K (29) 


This resembles the formula for the magnetic drift velocity in the 
presence of a force field - VK , but is actually not the same thing 
(except in the limit J = 0 , when the particles are confined to the equa- 
torial plane) since V is an averaged drift velocity, averaged over the 
entire field line and then projected onto the equatorial plane. 
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Equation (29) provides a method for calculating V ; for instance, in 
the case of the dipole field or in other axisymmetric fields where ^ = af 
one obtains, by resolving B into x 7 ^ , the relation 

(i+ 1.1254 x - 0.025664 X^ 5 +- 0.0727 + 2 X ^/5 I~ ) (30) 

eaL v 

which also follows from (25)® This can be evaluated by means of (26) and 
provides a quick derivation of the averaged magnetic drift velocity 
(electric drifts can be simply added to it) which in general requires a 
lengthy calculation. 

If the values of (j, |u) are set, K is a function of ( ol, ) and 
can be chosen to serve as an Euler potential. Let then u(ot,^ ) be the 
conjugate Euler potential, i.e. let it satisfy 

B = V K x 7 u( ) ( 51 ) 

The mean rate at which u changes can be evaluated in the equatorial 
plane and is, in the absence of explicit time dependence 


<TX> - V - Vu = e/c = constant (52 ) 

Taken together, K and u describe the mean motion of particles, mapped 
into the equatorial plane, in a simple intuitive way. By (29) 


V-VK = 0 ( 53 ) 

i.e. the particle stays on a constant value of K (this no longer 
holds if K contains direct dependence on t), so that lines of constant 
K in the equatorial plane resemble rays along which particles are guided. 

On the other hand, by (52) lines of constant u play the role of "wavefronts ", 
since a group of particles starting initially from a given value of u 
will at al l times share the same value of u , in the particular represen- 
tation of u which was initially chosen (the choice of u is not unique. 
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and a group of particles lined up along a line of constant u in one 
representation might no longer do so when another one is used)* 


These "pacing functions" (they pace the rate at which particles advance 
along lines of constant K) will be discussed in a separate paper. The 
"conjugate potential" of a given electric potential 4 * ( f Stern , 
1974 a J represent the special case of u when W = 0 and ^/^t = 0 
in equation ( 7 ) • 


According to Northrop and Teller, the "bounce time" T for one 
complete oscillation between mirror points satisfies 

T _1 = ^K/T>J 

= ( ■BK/'DX) X/T> J) (34) 


By ( 15 ) X depends on J through l(y) (the denominator there 
is related to the first invariant p) and explicitely one then finds 


X = j/2mvaly 


(35) 


(this may be compared to eq. 9 > where I has attained its maximal 
value Iq )• Thus 

T" 1 = ( p|g^|/ 8 m a 2 L 5 ) 2 ( 1.3505 - 0.040567 X 1 / 5 + 

4- 0.16777 X 2 ^ + 2 X / 3® ) (36) 

Note that the averaged equatorial drift velocity V and bounce period 
T given by equations (30) and ( 36 ) respectively reduce to the 
expressions derived by Schulz [Mi] in the limits of J ~ 0 and |jl = 0 
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CAPTIONS TO FIGURE AND TABLE 


Figure 1 


A plot of X (defined in equation 15 ) against the 
equatorial pitch angle > . 


Table 1 - Values of the product y 2 f (x) (given in equation 20) 

using the approximation ( 19 ) for f(x), as a function of 
mirror colatitude 0^ and of y 2 . 



w 

Table 1 


9 m 

degrees 

> 

degrees 

y* 

y 2 f(x) 

5 

0.03 


1.0290 

10 

0.21 


1.0356 

15 

0.71 

0.0002 

1.0297 

20 

1.66 

0.0008 

1.0179 

25 

3«17 

0.0031 

1.0063 

30 

5.34 

0.0087 

0.9989 

35 

8.23 

0.0205 

0.9966 

40 

11.89 

0.0424 

0.9986 

45 

16.33 

0.0791 

1.0032 

50 

21.56 

0.1350 

1.0075 

55 

27.58 

0.2143 

I.0099 

60 

34.38 

0.3189 

1.0102 

65 

41.97 

0.4473 

1.0081 

70 

50.32 

0.5924 

1.0049 

74 

57.54 

0.7120 

1.0026 

78 

65.20 

0.821*2 

1.0010 

82 

73.23 

0.9167 

1.0001 

86 

81.54 

0.9782 

O.9998 
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